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• • . Abstract 

^ ■ Renormalization-group methods in soft-collinear effective theory are used to perform the 

^ , resummation of large perturbative logarithms for deep-inelastic scattering in the thresh- 

■ " " ' old region x — > 1. The factorization theorem for the structure function F2{x,Q'^) for 

X ^ 1 is rederived in the effective theory, whereby contributions from the hard scale 
and the jet scale Q^(l — x) are encoded in Wilson coefficients of effective-theory oper- 
ators. Resummation is achieved by solving the evolution equations for these operators. 
Simple analytic results for the resummed expressions are obtained directly in momentum 
space, and are free of the Landau-pole singularities inherent to the traditional moment- 
space results. We show analytically that the two methods are nonetheless equivalent 
order by order in the perturbative expansion, and perform a numerical comparison up to 
next-to-next-to-leading order in renormalization-group improved perturbation theory. 



1 Introduction 



It is well known that fixed-order perturbation theory is not reliable for quantities involving 
several disparate scales. In such cases, higher-order corrections are enhanced by large loga- 
rithms of scale ratios. The standard solution to this problem is to split the calculation into a 
series of single-scale problems by successively integrating out the physics associated with the 
largest remaining scale. Perturbative logarithms are then resummed by renormalization-group 
(RG) evolution from the larger scales to the smaller ones. For collider processes, resummation 
is traditionally performed by other means, since it was not always clear how to systematically 
integrate out the physics associated with high scales in such cases. 

The simplest example of a high-energy process with a scale hierarchy which necessitates 
resummation is deep-inelastic scattering (DIS) in the threshold region. As the Bjorken scaling 
variable x — >■ 1, the invariant mass of the hadronic system produced in the decay, Mx = 
Q\/^—^ (neglecting the nucleon mass) 



becomes much smaller than the momentum transfer 
Q. The presence of the two scales is manifest in the QCD factorization theorem P3I21|3J 



Fr(x,g2) = i^(Q^/i)g2 f'-j(Q 

J X Z ^ 
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for the non-singlet part of the structure function F2{x,Q'^). The result is valid in the 
threshold region at leading power in Mx/Q"^ ^ (1 — x) and Aqqj^/Mx- As long as Mx 3> 
Aqcd, both the jet function J{Mx,fj,) and the hard function H{Q'^,n) can be evaluated in 
perturbation theory, whereas the parton distribution function is a non-perturbative 

object. The result for the hard function involves single and (Sudakov) double logarithms of the 
form a" ln'"((5//x), with m < 2n, while the integral over the jet function produces logarithms 
a" ln™'(Mx//i). Irrespective of the value of the renormalization scale /i, the fixed-order result 
contains large logarithms. 

Traditionally, the resummation of these logarithms is performed in moment space. The 
threshold region of small Mx is probed by large-A^ moments. The relevant scale in Mellin 
space is Q/\fN, so that the large perturbative logarithms depend on the moment parameter 
A^. In pi, IE| it was shown that these logarithms can be absorbed into a resummation exponent 
Gat, defined by integrals over two radiation functions Aq{as 



and Bq{as 



Gn{Q = / dz 
Jo 



— Agiasik)) + Bg [asiQVT^ 



(2) 



The functions Ag and Bg are determined by matching with results from fixed-order perturba- 
tion theory and are currently known at three-loop order, enabling a nearly complete threshold 
resummation to next-to-next-to- next-to-leading logarithmic (N^LL) accuracy 4j. The re- 
summed momentum-space structure function F2{x,Q'^) is obtained from the moment-space 
expression by an inverse Mellin transformation. 

This approach to threshold resummation has several drawbacks. The first is related to 
integrations over the Landau pole in the running coupling. These occur twice: once in the 
integrals over the functions Ag and Bg in the resummation exponent, and once again when 
the inverse Mellin transform is taken to obtain results in momentum space. To perform the 
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resummation one needs to specify a prescription for how to deal with these poles. Various 
methods have been proposed in the literature, such as the "minimal prescription" or the 
"tower expansion" [Uj. The difference between these prescriptions is a power-suppressed effect. 
Since factorization theorems do receive power corrections, this does not appear as a problem 
at first sight. However, as discussed in [7j, the Landau-pole singularity in the resummed 
expression can induce large unphysical power corrections. In the example of the Drell-Yan 
process, the ambiguity in the threshold resummation amounts to a power correction of order 
Aqcd/^x, while the physical power corrections to the process scale as Aqqj-,/M|-. The fact 
that resummations with RG methods [3 |HI do not involve Landau-pole ambiguities illustrates 
that these effects do not have a direct physical interpretation. In particular, a Landau-pole 
ambiguity does not necessarily imply the presence of a commensurate renormalon ambiguity 
[Zj. Further drawbacks are that in the traditional resummation formalism the separation of 
contributions from the hard and jet scales is not transparent, and while the function Ag has a 
clear interpretation as the cusp anomalous dimension familiar from the renormalization theory 
of Wilson lines PlCni; the function Bg is not easily identified with a field-theoretical object. 

In this paper we use RG techniques to perform the resummation of perturbative loga- 
rithms directly in momentum space. The starting point is the factorization formula ([Q), which 
we rederive using soft-collinear effective theory (SCET) jTT| ^[ E[ ■ In this framework, 
the hard function H and the jet function J are matching coefficients. The hard function 
arises from a first matching step, in which the electroweak current is matched onto a corre- 
sponding effective-theory current operator. In a second step, the partons associated with the 
hadronic final state are integrated out, giving rise to the jet function. Threshold logarithms 
are resummed by solving the RG equations for these matching functions, using techniques 
presented in ^]. Existing results from higher-order perturbative calculations enable us to 
perform the matching and resummation up to next-to-next-to-leading order (NNLO) in RG- 
improved perturbation theory, corresponding to the N^LL approximation in the standard 
approach. We show that the results obtained in momentum space are formally equivalent 
to the more familiar moment-space results by deriving a formula which connects them order 
by order in perturbation theory. However, integrals over the Landau pole never appear in 
our momentum-space formulation, and the effective-theory matching functions and anoma- 
lous dimensions have a clear field-theoretical interpretation. Furthermore, we obtain a simple 
analytic expression for the resummed structure function, while the Mellin inversion which is 
necessary in the traditional approach can only be performed numerically. As a result, it is 
straightforward to match our resummed expressions onto fixed-order calculations valid outside 
the threshold region. Finally, we stress that our approach to resummation in x-space is free of 
the pathologies related to large unphysical power ambiguities found in In fact, it exhibits 
a better apparent perturbative convergence than the conventional approach. 

In the context of SCET, the generic factorization formula for DIS has been discussed 
previously in [16j, while the case a; — > 1 has been studied in [T7| UHl CHI 1201 1211 122] ■ These 
papers make conflicting statements about the factorization properties of DIS in the endpoint 
region. Most of the differences are resolved after observing that, near the endpoint, the 
parton distribution function receives contributions from two distinct non-perturbative modes. 
While the two modes cannot be factorized perturbatively, their presence must nonetheless be 
taken into account to correctly translate the effective-theory result into the QCD factorization 
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theorem (jT)). In [7| [T7| UHl HH I2S| , the resummation was performed by solving the RG equations 
for the moments. This avoids Landau-pole singularities in the exponent, but as we show here, 
it is possible to solve the equations directly in momentum space. 

While threshold resummation in DIS is of limited phenomenological importance, it is a 
relatively simple process for which the perturbative results needed in our calculation are known 
at NNLO. For this reason, it provides an especially instructive example with which to develop 
our resummation formalism. However, many other processes fulfill factorization theorems of 
the same structure, in which the rate factorizes into a hard contribution times a jet function 
convoluted with a nonperturbative matrix element, and our formalism also applies to these 
cases. An example is heavy-particle production near threshold, which includes the Drell- 
Yan process in the limit where the invariant mass of the produced lepton pair is close to 
the center-of-mass energy in the collision, as well as Higgs production in the same kinematic 
region. Another example is provided by inclusive S-meson decays in the endpoint region. Our 
final result for the resummed DIS structure function is very similar to the factorized expression 
for radiative B Xg'j decay as derived in |211I2S1- In fact, both processes involve the same 
jet function, given by the quark propagator in light-cone gauge. 

The outline of the paper is as follows. In Section |21 we use SCET to obtain the QCD factor- 
ization formula for DIS near the endpoint, providing a translation between the effective theory 
and standard discussions. In Section El we work out the technique for threshold resummation 
in momentum space, derive a compact expression for the factorized structure function, and 
give results for the perturbative matching coefficients valid to NNLO in perturbation theory. 
In Section |3] we convert our results to moment space and show how they are connected to 
those obtained in the standard approach. The Appendix gives the perturbative expansions of 
the RG functions used in our analysis. 



2 Factorization in DIS 

In this section we derive the QCD factorization formula for the non-singlet DIS structure 
function F^%x,Q^), using the technology of SCET [m [H HSJ HH. We consider DIS of 
electrons off a nuclear target, e~ + N{p) e~ + X{P), as illustrated in Figure ^ All non- 
trivial hadronic physics is encoded in the hadronic tensor 

W>"'ip,q) = i f c/^xe'^-" {N{p)\T{f>'{x) ^(0)} \N{p)) 



averaged over the nucleon spin. Here J'^ = ip'j'^ip is the electromagnetic current. The scalar 
functions Wi can be expressed in terms of the kinematic invariants 

Q^ = -g^ ^ = ^, (4) 

2p ■ q 

where q = P — p is the momentum of the virtual photon, and x is the Bjorken scaling 
variable. For simplicity we focus on the flavor non-singlet component of the cross section. 



3 



Figure 1: Kinematics of DIS. 



which is insensitive to the gluon distribution in the nucleon. It can be obtained by taking the 
difference of the DIS cross sections for scattering off different target nuclei. 

The first step in analyzing the hadronic tensor is to identify which momentum regions give 
non-vanishing contributions to the Feynman diagrams for W'^'^ in perturbative QCD. Each 
region is represented by a set of fields in SCET. The identification of regions can be done 
in any Lorentz frame. The number of modes and their relative scaling is Lorentz invariant.^ 
Two particularly convenient reference frames are the target rest frame, where = (m, 0, 0, 0) 
with m the nucleon mass, and the Breit frame, where the virtual photon carries momentum 

= (0,0,0,(5). The two frames are related to each other by a Lorentz boost along the 
z-direction. Introducing the light-cone decomposition 

p'' = ■ p) Y + ■ p) Y + p± = p+ + p'^ + p± , (5) 

where = (1,0,0,1) and n'^ = (1,0,0,-1) are two light-like basis vectors (n ■ = 2), a 
generic momentum p'^^^^ = p'^ + pt + p'^ in the target rest frame transforms into Psreit — 
p^ + pt +Px Breit frame, with the rapidity t] of the boost given by 



We shall discuss the different regions in the Breit frame, where the final-state hadronic jet 
moves along the 2;-direction, while the target nucleon moves in the opposite direction. The 
light-cone projections of the relevant momenta are (all perpendicular components vanish by 
choice of coordinates) 

n - q = ~Q , n - q = Q , 

Q rri^x _ _„ m^x 

n ■ p = me' = 1 — — — h..., n ■ p = me ' = — — h . . . , 

X Q Q 

^ ^ 1 — X m^x „ ^ m'^x 

n-P = Q —— + -^ + ..., n-P = Q + -^ + ... , (7) 



^We disagree with the claim of jl7) that fewer momentum regions contribute in the target rest frame than 
in the Breit frame. 
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where the neglected terms are of order m'^x^/Q'^. The invariant mass Mx of the final-state 
hadronic jet is given by 

X X 

In the last step we have used that for an inclusive process the jet mass must be much larger than 
m ~ Aqcd- Otherwise, the cross section cannot be analyzed using short-distance methods. 

While the momentum of the virtual photon is fixed by kinematics, the final-state jet and 
target nucleon consist of jets of near on-shell partons, whose momenta have scalings consistent 
with the relations above. We introduce a small expansion parameter A ~ m/Q ~ Aqcd/Q 
and quote the components {p+,p-,p±) of parton momenta in units of Q. Assume first that 
X = 0{1) is not very close to 1. Then, for the purposes of power counting, it follows that 
g ~ Q (1, 1, 0), p ~ Q (1, A^, 0), and P ~ Q (1,1,0). The partons making up the initial and 
final hadronic states have generic scalings 

target nucleon: pc ~ Q (1, A^, A) (anti-coUinear) , 
final-state jet: Ph ^ Q (1, 1, 1) (hard) , (9) 

where the term "anti-collinear" refers to collinear fields propagating in the negative ^;-direction. 
These relations change for the special case where x is close to 1, such that e = 1 — x becomes 
parametrically small. The momentum of the final-state jet now scales like P ~ Q{e,l,0). 
While the valence quark in the target nucleon struck by the photon still carries an anti- 
collinear momentum Pvaience ~ Q(l5A^,A), the remaining partons in the nucleon now have 
total momentum scaling like p — Pvaience ~ Q (e, A^, A). Consequently, the partons making up 
the final-state jet and the target remnant jet have momenta scaling like 

final-state jet: phc ~ Q (e, 1, v^) (hard-collinear) , 
target remnants: Psc ~ Q (e, A^, v^A) (soft-coUinear) . (10) 

In these relations, the scaling of the perpendicular momentum components follows from the 
requirement that the individual partons be nearly on-shell. The terminology for the "hard- 
collinear" and "soft-collinear" modes follows |26j and [2ZI. In the traditional literature on 
factorization in DIS ^1 |2l El the soft-collinear modes were referred to as "soft" . Relation (jHJ 
implies that A^ <^ e ^ 1, and there is no need to specify the relative scaling between e and A in 
more detail. All that matters for the factorization analysis is that p1^ ~ Q^e is a perturbative 
scale, while ~ Q^X^ is not. 

The discussion of factorization for the generic case, where x = 0{1) but not very close 
to 1, is straightforward 16j. Hard modes are described by QCD, whereas the anti-collinear 
partons making up the target nucleon can be described in SCET. There is no need to include 
any other modes, since the only relevant regions are hard and anti-collinear. In interactions 
of the anti-collinear fields with hard fields, only the large plus components Pc+ ~ Q of the 
anti-collinear momenta should be kept at leading order in power counting. Correspondingly, 
the anti-collinear fields must be multipole expanded about a;_ = (n ■ x)n/2, i.e., (psi^) = 
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4>c{x-) + . . .. Integrating out the hard modes by matching onto SCET yields an expression 
for the discontinuity of the hadronic tensor of the form 

-lmW>^^= ^^C(Q^x/e,/x) / ^e-^«-^*(iV(p)I^M[tn,0]7^^7^^(0)|iV(p)), (11) 

TT Jx X J ZTT Z 

where C = 5(1 — x/C,) + 0{as) is a matching coefficient in the effective theory, and the object 
[tn, 0] is a straight Wilson line along the n light-cone. We have used that the SCET Lagrangian 
for a single collinear sector is equivalent to the QCD Lagrangian jT2I in order to replace the 
SCET fields by the usual QCD fields. The identification of the nucleon matrix element with 
the QCD parton distribution function is then automatic (see relation ()23p below), and one 
arrives at the standard factorization formula. 

The derivation of the factorization formula for x — > 1 is more complicated. It involves a two- 
step matching procedure similar to that used for inclusive semi-leptonic and radiative B decays 
in the endpoint region [211 EHl EEl EB] • In a first matching step, hard modes are integrated out 
by matching QCD onto a version of SCET containing hard-collinear, anti-collinear, and soft- 
coUinear fields. We will refer to this intermediate effective theory as SCET(/ic, c, sc) for short. 
The matching function associated with this first step is the hard coefficient Cy. Because the 
sum of a hard-collinear momentum and an anti-collinear momentum has an invariant mass 
(Phc + Pc)'^ ~ and must be counted as hard, the intermediate effective Lagrangian does not 
contain vertices coupling the hard-collinear fields to anti-collinear ones. These fields interact 
only through the exchange of soft-collinear "messenger" fields. However, the soft-collinear 
modes can be decoupled from the hard-collinear ones by means of a field redefinition. After 
this decoupling, it is possible to integrate out the hard-collinear scale by matching onto a 
low-energy theory SCET(c, sc) involving only anti-collinear and soft-collinear modes. The 
matching function associated with this step is the jet function J. Having integrated out the 
perturbative modes, the final step is to evaluate the matrix element of the remaining operator 
defined in the low-energy effective theory. An important part of the factorization analysis is to 
show that this matrix element is equivalent to the QCD parton distribution function evaluated 
in the limit studied e.g. in Ej. We will show that this is indeed the case, and 

that the soft-collinear modes play an important role in this identification. 

The appropriate Lagrangian for SCET(/ic, c, sc) is a generalization of the effective La- 
grangian for collinear and soft-collinear fields derived in |27| 1^ . It contains hard-collinear 
quark and gluon fields ^hc and A^c, anti-collinear quark and gluon fields and A^, and 
soft-collinear quark and gluon fields O^c and Ag^- The hard-collinear fields move along the 
2;-direction, and hence yi^hc = 0. The anti-collinear and soft-collinear fields move in the 
opposite direction, so = and ^9sc = 0. The two collinear sectors can only interact 
via soft-collinear exchange, and at leading power only soft-collinear gluons are involved in 
these interactions. The corresponding effective Lagrangian at leading order in the expansion 
parameters e and A is j2Zl 1201 

— ^ — 1 

'^^scet(2/) = ^hc - [in ■ Dhc + gn ■ Asdv-)] ^hc - ^hc iphc^ - — — Tr~ ^Phc± ^hc 

— % — IJL \ 

-I- - [in ■ Ds + gn- Asc{y+)] ^5 - ^cipc± —Jj- ^PcL 
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+ pure glue terms + soft-collinear Lagrangian , 



(12) 



where all fields without position argument are to be evaluated at the point y. The effective 
Lagrangian is invariant under a set of hard-coUinear, anti-coUinear, and soft-collinear gauge 
transformations, whose precise form can be found in p7| IHOj. 

An important property of the SCET Lagrangian is that soft-collinear gluons can be de- 
coupled from the hard-coUinear and anti-collinear fields through field redefinitions involving 
Wilson lines |121 ITfj . This decoupling is essential for the factorization analysis below. Di- 
agrammatic factorization proofs also rely on the decoupling of "soft" gluons from coUinear 
fields. The underlying physics is that soft gluons couple to collinear partons through eikonal 
vertices, a feature explicit in the SCET Lagrangian (jl2j) . 

2.1 Matching of the current 

The first step in the factorization procedure is to integrate out hard fluctuations by match- 
ing QCD onto the intermediate effective theory SCET(/ic, c, sc). The kinematic restrictions 
implied by the limit x ^ 1 simplify this first matching step. Since we are dealing with the 
region of phase space where the final-state jet is hard-collinear, there are no contributions to 
the hadronic tensor where the anti-collinear partons at points and x are connected by hard 
gluons. It is therefore sufficient to integrate out hard fluctuations at the level of the elec- 
tromagnetic current. Time-ordered products of two currents are not needed until the second 
step. 

We match the QCD current J^{x) = {'ijj'j'^ijj){x) onto a current in SCET containing a 
hard-collinear quark and an anti-collinear anti-quark. The form of the resulting operator is 
dictated by gauge invariance. The appropriate matching relations for the QCD fields are 

i^hcix) ^ iwUhc)ix) , ^ {W^^,){x.) , (13) 

where Whc is the hard-collinear Wilson line 

Whc{x) = P exp (^ig J dsn- A^^x + sn)^ (14) 

along the n-direction, and Wc is the analogous anti-collinear Wilson line along the n-direction. 
The multipole expansion in (fT^ requires some explanation. In the hadronic tensor Q, the 
points and x are connected by a hard-collinear jet propagating through a cloud of soft- 
collinear partons. This implies that the position argument x scales as a hard-collinear quantity, 
X ~ (l,e~^,e~2). It follows that not all components of the anti-collinear and soft-collinear 
momenta must be kept in the calculation of Feynman graphs in the effective theory. The minus 
and perpendicular components of anti-collinear and soft-collinear momenta are much smaller 
than the corresponding components of hard-collinear momenta and so should be expanded 
out. On the other hand, the large plus component n • p ~ Q of the target nucleon is canceled 
by the momentum component n ■ q oi the current and turned into a momentum component 
of order eQ, which is of the same order as the plus component of a hard-collinear or soft- 
collinear momentum. For this reason, it would be incorrect to set x_ = in the argument of 
the soft-collinear fields entering the effective current operator, even though this is the correct 
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multipole expansion of Lagrangian interactions between soft-collinear and anti-collinear fields 
|27j . Therefore, when matching the current operator onto SCET, one must multipole expand 
both the anti-collinear and soft-collinear fields about 

The two expressions in (fT^ are invariant under hard-collinear and anti-collinear gauge 
transformations, while under a soft-collinear gauge transformation both composite fields, 
Wl^^hc and W^^c, transform into Usc{x^) times themselves. Thus, at tree level the gauge- 
invariant matching relation for the current is 

^ {^,W,){x.) ll{wUhc){x) . (15) 

Only a single Dirac structure is possible for massless quarks. Beyond tree level the matching 
relation at leading power gets generalized to (see the analogous discussions in 1^) 

(^7^^)(x) ^ I dtCvit,n-q,fx) (^,W,) (x ^) ^^iwU,,) (x + tn) 

= Cvi-n-qn-P,fx){^,W,)ix^)^^^iwUhc)ix). (16) 

In the first line we have used that n ■ d derivatives of hard-collinear fields are unsuppressed 
in SCET power counting, allowing for arbitrary displacements of these fields along the n 
light-cone. In the second line, the object P is the hard-collinear momentum operator, and 
the Wilson coefficient Cy is the Fourier transform of the position-space Wilson coefficient Cy 
appearing in the first line. In the case at hand, the relevant components —n ■ q k, n ■ P Q 
are fixed by kinematics (see the relations ((Tj)), and so we may write CyiQ'^^lj) for simplicity. 

2.2 Matching of the hadronic tensor 

The next step in the matching procedure is to evaluate the hadronic tensor in the intermediate 
effective theory. Inserting the SCET current (pi)|l into (jH)), we find the leading-power expression 

iy^^(p,g) ^ \Cv{Q^^^)?^ J d^xe"^'^ (17) 

The interactions of soft-collinear gluons with hard-collinear fields in (jl2j) can be removed by 
the field redefinitions (HI IZZ] 

ac(x) - S^{x^) d?(a:) , Alix) - 5„(x_) A',f\x) SHx.) , (18) 

which imply (WUhcKx) ^ 5„(a;_) (ivif efi)(a;). Here 

Sn{x) = P exp (^g J dsn ■ Asc{x + sn)^ (19) 

is a soft-collinear Wilson line along the n-direction. The redefined hard-collinear fields with 
superscripts "(0)" are decoupled from soft-collinear fields and thus interact only among them- 
selves. After the field redefinition the hadronic matrix element in ()17|) factorizes into a vacuum 
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(a) (b) (c) 

Figure 2: Examples of diagrams involving anti-collinear gluon exchange. The dashed (dot- 
ted) lines represent anti-collinear (hard-collinear) quark lines. The wavy lines represent the 
electromagnetic currents. In graph (a) the anti-collinear gluon is part of the initial-state (nu- 
cleon) jet and the final-state propagator is hard-collinear, as required in the effective theory. 
In graphs (b) and (c) the anti-collinear gluon is part of the final-state jet, whose invariant 
mass then becomes hard. Graphs (b) and (c) are therefore not part of the effective-theory 
representation of the hadronic tensor as a; — > 1. 



matrix element of hard-collinear fields and a nucleon matrix element of anti-collinear and soft- 
collinear fields. In the second matching step, we "integrate out" the hard-collinear fields, which 
can be done using perturbation theory because the hard-collinear scale is a short-distance scale, 
(5^(1 — x) ^ ^QCD- Since in a single (hard-)collinear sector SCET is equivalent to full 



QCD the vacuum matrix element of hard-collinear fields can be rewritten in terms of the 
QCD matrix element ^Slj 



|0) 



(20) 



The object W{x) denotes a Wilson line analogous to (fT^ but with gauge fields in full QCD. 
Color indices are suppressed; the correlator is proportional to the unit matrix in color space. 
We define the jet function through the imaginary part of J as (see e.g. j32) 



J(p^, /i) = — Im 

TT 



^^(p^ /u) 



(21) 



The jet function has support for > 0. 

At this point it is important to emphasize a subtlety related to the matching of the forward- 
scattering amplitude in full QCD onto operator matrix elements in SCET. The anti-collinear 
composite fields W^^c and ^cWc in (fTTjl are not allowed to communicate via anti-collinear 
particle exchanges, but only through exchanges of soft-coUinear partons. The exchange of anti- 
collinear particles between the two currents in (jHl) is kinematically forbidden in the region x 
1, as this would lead to a final-state invariant hadronic mass Mx ~ Q- Since the intermediate 
state is hard instead of hard-collinear, diagrams such as those shown in Figure Efb) and 
(c) are not part of the effective-theory representation of the hadronic tensor in the region 
X — > 1. Such "forbidden" graphs are nonetheless generated (and yield non- vanishing results) 
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if the SCET Feynman rules used for the matching of the electromagnetic current are naively 
applied to the hadronic tensor. We can construct a set of Feynman rules appropriate for the 
hadronic tensor by introducing different anti-collinear fields for the "in" and "out" states in 
the forward-scattering amplitude, and restricting interactions between the two anti-collinear 
sectors to soft-coUinear exchange. These effective-theory Feynman rules produce graphs such 
as that in Figure Efa), but not those in Figure |2fb) and (c). For simplicity of notation, we 
will suppress the "in" and "out" labels on the anti-collinear fields, but one must make this 
distinction when evaluating the hadronic tensor in the effective theory. 

After integrating out the hard-coUinear fields, the resulting nucleon matrix element in the 
low-energy effective theory can be reduced to 

{N{p)\ ^n(x-) 1^ \ 7i ^i(o) {wlim 

= -{Nip)\ i^mi^-) [x-,0]sc{9T - ^6^75) f (W^^m \Nip)) , (22) 

where 0]sc = Sn{x_) 5*^^(0) is a straight Wilson line of soft-coUinear gluon fields along the 
n light-cone. In the second line we have defined the objects g^^ = g^'^ — \ {n^n'^ + n^rf) 
and e^^ = \t^^"^na_np, and also used that = 0. The anti-symmetric structure vanishes 
after averaging over the nucleon spin. The appearance of the symmetric structure g^ implies 
the Callan-Gross relation W2 = Wi at leading power and to all orders in perturbation 
theory. Hereafter, we thus focus on the structure function Wi. 

Consider now the standard definition of the quark distribution function in QCD [33j . 

^ ZTT J~od I 

where [tn, 0] is a straight Wilson line of gauge fields in full QCD, and the superscript "ns" 
indicates the flavor non-singlet component of the distribution function. In the Breit frame, 
where the proton moves along the n-direction, ■?/' and can be considered anti-collinear fields. 
For generic values of ^ these fields carry only a portion of the proton's longitudinal momentum. 
The remaining portion (1 — ,^) n ■ p is still large and can be shared between other anti-collinear 
partons exchanged between the two points and tn. A different picture is called for in the limit 
^ — >■ 1, where the anti-collinear valence quarks i\} and i\) carry almost all of the longitudinal 
momentum [H]. In this case, the residual momentum component (1 — ^)n ■ p is small, and 
the remaining partons are soft-collinear. Each valence quark is described by an anti-collinear 
jet propagating through the soft-collinear cloud made up of the remaining partons. The two 
anti-collinear jets communicate through soft-collinear gluon exchange only. 

The distinct roles played by the valence and remaining partons as ^ — 1 make it appro- 
priate to introduce an effective field-theory description for the parton distribution function, in 
which it is matched onto an operator involving anti-collinear and soft-collinear fields in SCET. 
The most general, gauge-invariant form the relation (|23|) can be matched onto in the ^ ^ 1 
limit reads 

= ^ jyte~^^'-^{N{p)\ {^,W,){tn) [tn, 0]., | (l^t^,)(0) \N{p)) . (24) 
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It is understood that the anti-coHinear fields located at the points and tn interact only via 
soft-collinear gluon exchange. Both (|24j) and the QCD matrix element (|23p depend on the 
single invariant = m?, so there is no non-trivial hard matching coefficient. The matrix 
element is precisely the object we encountered in ((221) • We can use this correspondence 
along with some simple algebra to find 

Wi = \Cv{Q\^^)\^^ I d{n-k)n-qj{q^ + n-kn-q,ii)(PlH'^,ii) . (25) 
J ^ ■ p 

The structure function F^^^x^Q"^) equals e^x^ImWi, where the Cq are quark electric 
charges. Inserting the definition of the jet function (j2^, and recalling that q^ = —Q^ and 
n ■ pn ■ q = Q'^/x + power corrections, we obtain the final result for the factorization formula 

Fn^, Q') = E \CviQ', /x) r Q' tdi J{Q' /x) ^^) . (26) 

q J X ^ X ' 

This formula is valid to all orders in perturbation theory and at leading power in (1 — x) and 
^qcd/^x- '^^^ argument of the jet function takes values between and M|^, where the total 
jet invariant mass was given in (jSJ. The equivalent form (0) is obtained by substituting ^ = 
x/ z. At tree-level, this formula evaluates to the familiar parton-model expression F2^{x, Q"^) = 
Eqelx<P^%x). 

Relation ()26|) is the standard form of the QCD factorization formula for the DIS structure 
function in the limit x — 1 |21 EI , which we have derived here using SCET. We hope our 
derivation helps resolve some of the disagreements in the literature. Soft-collinear messenger 
modes obviously play a crucial role in the derivation, as the parton distribution function 
at large ^ is defined in terms of these fields. The proper effective-theory description of the 
parton distribution function thus requires two distinct non-perturbative modes. This element 
is missing from fri\ 120] , where it was argued that only one non-perturbative mode is needed, 
either because the soft graphs vanish in the Breit frame calculation, or because the effective- 
theory formulation in the target rest frame involves only one non-perturbative mode from the 
beginning. Although we disagree with these statements (the second of which would violate 
reparameterization invariance in the effective theory), our explicit one- loop results agree with 
those derived in these papers. This is because our findings imply that parton evolution in 
the endpoint region can be described simply by taking the x — 1 limit of the Altarelli-Parisi 
splitting functions, which is effectively what was done in the calculations of jT^j . Our explicit 
one-loop results also agree with those in |18j, where the power counting e = l — x~A = 
^qcd/Q was adopted. While this counting is possible and natural in view of the hierarchy 

<^ e ^ 1, it does not imply that the soft-collinear scale m^(l — x) depends on the scale 
Q, and the presence of this scale does not translate into non-perturbative Q-dependence 
in the parton distribution function, as was suggested in Finally, we have shown that 

the soft-collinear contributions are precisely such that they can be absorbed into the parton 
distribution function. We therefore do not confirm the claims of soft contributions outside the 
parton distribution function made in jTH] . The same conclusion as ours was reached in [211 122] , 
where it was argued that the infrared divergences due to collinear and soft emissions can be 
absorbed in the standard QCD parton distribution function, although |21] did not discuss 
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tn 








Figure 3: Soft-coUinear Wilson line Wc- 

how to obtain this result in the effective theory. In j221 it was claimed that there is a double- 
counting problem in SCET, which must be remedied by making certain soft subtractions from 
the coUinear matrix element. We have shown here that there is no such problem. Double 
counting occurs only if one fails to notice that collinear emissions such as those shown in 
Figure Efb) and (c) must not be included in the effective-theory calculation near the endpoint. 
Similar to ^\ , the discussion in [221 fails to distinguish the virtualities of hard-coUinear and 
anti-collinear modes, and it overlooks the fact that the smallest scale in the problem is not 
Q^{l-xf but m'^{l-x). 

In the next subsection, we will emphasize the importance of soft-coUinear Wilson loops in 
determining the RG properties of the effective theory. 

2.3 Decoupling transformation and cusp singularities 

An important step in the derivation of the factorization formula ()26|) was the identification of 
the parton distribution function for ^ — > 1 with the SCET matrix element on the right-hand 
side of fl24j) . We can simplify this relation further by decoupling the soft-collinear gluons from 
the anti-collinear fields with the help of the field redefinitions 

Uy) - Sn{y+) #(2/) , A^M - Sn{y+) Af\y) Sl{y+) , (27) 

where the soft-collinear Wilson line Sn is defined in analogy with but with n replaced 
by n. Above, j/ is a generic argument of a term in the SCET Lagrangian. The multipole 
expansion of the soft-collinear fields about must be done everywhere except at the location 
of the current, where x scales as a hard-collinear (not anti-collinear) quantity, see above. At 
this one point, we have instead 

ec-(x-) -5^(.T_)#(a;_), A^^ix.) ^ Sn{x.) Af\x.) Slix.) . (28) 

It follows that 

= V- r^^e-^*"-"(Ar(p)| (e?%f )(tn) ^H^c(t) (W^F^f )(0) \Nip)) , (29) 

where 

Wcit) = (0| Slitn) [tn,OUSn{0) |0) = (0| SUtn) S„itn) 31(0) 3^(0) |0) (30) 
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describes a closed Wilson loop consisting of the junction of a Wilson line extending from 
— oo to along the n-direction, a finite- length segment from to tn along the n-direction, 
and another Wilson line from tn to — oo along the n-direction, see Figure El Anti-collinear 
virtual particles can be exchanged inside the brackets {wP^^i°^) and {^i^^W^) in ^ but not 
between them, see Figure El These exchanges give rise to non-perturbative renormalization 
factors Z 2(171, fi) u{p), where u{p) is an on-shell spinor, and the only invariant is = m? . The 
above formula then reduces to 

= ^K/^) ^ r dte-'^'--^Wc{t) . (31) 

This form of the parton distribution function for — > 1 coincides with eq. (2.12) of [S] (where 
our factor Z is called H). Performing the same decoupling transformation (j28j) on the SCET 
current (fTB|) yields 

(^7'^V^)(x) - CviQ^f^) (#iyf )(a;_) 4(^-) S„ix-) Y{wtTd^)i^) ■ (32) 

Once again, the soft-collinear fields reside in a closed Wilson loop, this time extending from 
— 00 to the point x_ along the n-direction, and returning to —00 along the n-direction. 

The appearance of soft-collinear Wilson loops in the SCET representation of the parton 
distribution function (j3ip and the electromagnetic current (|32p determines their RG properties. 
In both cases, ultraviolet singularities related to Sudakov double logarithms are governed by 
the so-called cusp anomalous dimension, Fcusp, which is a universal quantity of perturbative 
QCD Oinilini- We will confirm this structure in our exphcit calculations below. 

2.4 Power corrections 

Our focus in this paper is on the leading-order factorization formula (j^UI) . but it is important 
to keep in mind that this result receives power corrections in the small parameters e and A. A 
discussion of power corrections using SCET has been carried out for the closely related case 
of semi-leptonic inclusive B decay in the endpoint region, where it was found that the power 
corrections factorize order by order in l/rub P^ IH^ IHUj. It is evident that the same procedure 
applies here, so the effective theory offers a systematic tool to calculate power corrections for 
DIS near the endpoint. 

The procedure involves the two-step matching used in the leading-order case. A systematic 
treatment of the dominant power corrections requires to match the SCET Lagrangian and 
electromagnetic current at subleading order in e and A. We will not perform this complete 
matching here, but instead hmit ourselves to a qualitative discussion of the two types of power 
corrections: those to the jet function, which can be calculated perturbatively, and those to 
the parton distribution function, which cannot. We give examples of each, and explain how 
to obtain them with effective field-theory techniques. 

Corrections to the jet function first appear at order (1 — x) asifii) and can be calculated 
perturbatively. To identify the full set of such corrections, one must match the electromagnetic 
current onto SCET up to order e. Time-ordered products of the subleading SCET currents 
containing extra hard-collinear fields compared to the leading-order result (fT^ build up a 
set of power-suppressed jet functions convoluted with the leading-order parton distribution. 
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Figure 4: Examples of subleading time-ordered products in SCET which give rise to power 
corrections. Graph (a) leads to a perturbatively calculable power-suppressed jet function, 
while (b) and (c) lead to subleading parton distribution functions depending on three light- 
cone variables. 

As an example, consider the time-ordered product J'^(x) J'{0) involving two insertions of the 
0{y/e) suppressed current J' = n^^cWcWljphcL^hc- The relevant one-loop diagram is shown 
in Figure EJa). Decoupling the soft-collinear fields and factorizing them into the nucleon 
matrix element leaves the following vacuum matrix element of hard-collinear fields 

(0| T{(W^(?^rfie£))(a:) (^i?l!>i?)(0)} |0) . (33) 

The discontinuity of the Fourier-transformed matrix element defines a subleading jet function 
scaling as (1 — x) ^^(yUi). As in the case of the leading-order jet function, one could equally 
well calculate the matrix element in full QCD. This particular subleading jet function has 
been discussed previously in terms of a "non-local OPE" very close in spirit to SCET in [?7] . 
and was reconsidered in the context of SCET more recently in [TU] . 

Non-perturbative power corrections are given in terms of a basis of subleading parton dis- 
tribution functions. These are defined through nucleon matrix elements of power-suppressed 
SCET(c, sc) operators involving extra anti-collinear and soft-collinear fields compared to the 
leading-order matrix element (j^^ . To identify the complete basis, we would need to match 
both the current and the Lagrangians Cc+sc and Chc+sc to subleading order. An important 
simplification can be made by absorbing all time-ordered products involving insertions of the 
subleading Lagrangian terms Cc+sc into a redefinition of the leading-order parton distribution 
function 0(0- This can always be done, because such subleading parton distribution func- 
tions can only depend on a single light-cone variable ^ and are convoluted with the same jet 
function. On a more formal level, it amounts to not treating power-suppressed Cc+sc terms in 
the interaction picture, as is normally done. Identifying the remaining non-perturbative struc- 
ture then reduces to calculating power-suppressed time-ordered products involving subleading 
SCET currents or insertions of Chc+sc- 

The dominant subleading parton distribution function is related to the time-ordered prod- 
uct of the power-suppressed current arising from the multipole expansion of the anti-collinear 
fields, according to {WlQ{x) = [1 + ■ + i(a;x ■ d^f + . . .]{Wkc){x-). bmce x IS a 
hard-collinear quantity, x± ■ d±(j)c ~ (X/y/e) (pc- Thus the time-ordered product of J"\x) J(0), 
with 

J'\x) = \ (d ■ x^)\wUhc){x) (34) 
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results in a power correction of order A^/e. Because the power suppression comes from the 
current, the associated subleading parton distribution function depends only on a single light- 
cone variable, as was the case at leading order. 

Also interesting is the new multi-local structure related to interactions of soft-collinear 
fields with the jet. As an example, consider the time-ordered product 

j(°)t(a:)j(o)(0)/:Uc(l/)^Ue(^), (35) 

involving two insertions of the 0{\) suppressed Lagrangian term Cf^c+sc ~ Osc Ahc±C,hc- The 
relevant tree- level Feynman diagram is shown in Figure El^b). Performing the decoupling 
transformation and factorizing the soft-collinear fields leads to a tetra-local nucleon matrix 
element. An additional tetra-local correction of the same order is related to the emission 
of two transverse soft-collinear gluons, through the Lagrangian term C'l^c+sc ~ ihc^scL^hci 
as illustrated in Figure El^c). The parton distributions 0i(^i,^2,^3) defined by the scalar 
decomposition of the Fourier-transformed matrix elements appear in a three-fold convolution 
with a perturbatively calculable jet function, to be treated as power corrections of order 
A^Q;s(yu). For DIS, such multi-local hadronic structures have not yet been considered in the 
literature. 



3 Renormalization-group evolution and resummation 

The factorization formula for the DIS structure function derived in the previous section con- 
tains physics associated with different momentum scales factorized into a hard coefficient 
function Cy, a jet function J, and a non-perturbative parton distribution function 0°*^. These 
three objects depend on a scale /z at which the corresponding effective-theory operators are 
renormalized. The hard matching coefficient and the jet function need to be calculated using 
perturbative QCD. These calculations can be done at fixed order only when the scale is chosen 
so as to avoid large logarithms: the function Cy should be computed at a hard scale fih ^ Q, 
while the jet function should be computed at an intermediate scale Hi ~ Mx ~ QV^ ~ ^- The 
results of these calculations must then be evolved to the common scale ix in (j26j) by solving 
RG evolution equations. An advantage of the effective-theory formulation is that the RG 
equations can be solved directly in momentum space. The method has recently been outlined 
in [T3]; in this section we give a more detailed derivation, filling in the technical steps. 



3.1 Evolution of the hard function 



We begin by discussing the evolution of the hard matching coefficient Cy in |T6|) . At leading 
power there is a single gauge-invariant SCET operator the QCD current can match onto, and 
hence there is no operator mixing. The exact evolution equation takes the form 



d 



dlnfi 



rcusp('3^s 



In — + 7' 

/i2 



a. 



(36) 



The appearance of the cusp logarithm and its coefficient can be explained starting from 
using arguments presented in [21]. This term in the evolution equation is associated with Su 



dakov double logarithms. The remaining term, 7 , accounts for single-logarithmic evolution. 



15 



The exact solution to the evolution equation (|36p is 



2\ -"rC^/i./i) 



= exp 2S'(/i/i,/i) - a^v(/i/„/i) j Cv{Q'\iJ^h), (37) 

where /i/^ ~ Q is a hard matching scale at which the Wilson coefficient Cy is calculated 
using fixed-order perturbation theory. The Sudakov exponent S and the exponents a„ are the 
solutions to the differential equations 

S'(z/,/i) = -rcusp(as(/i)) In- , — ^ ar(i^,/x) = -rcusp(as(/i)) , (38) 



d m /i u a Inn 

and similarly for a^v, subject to the initial conditions S{v,v) = ar(z/, z/) = a^v(z/, z/) = 
at /i = z/. These equations can be integrated by writing d/dln/i = j3{as) d/das, where 
(3{as) = das/ din fi is the QCD /^-function. This yields the exact solutions OH] 

J /3(a) J (3{a') J /3(a) 

and similarly for the function a^v . The perturbative expansions of the anomalous dimensions 
and the resulting expressions for the evolution functions valid at NNLO in RG-improved 
perturbation theory are given in the Appendix. 



3.2 Evolution of the jet function 

The RG evolution of the jet function is more complicated. It was recently shown that the 
exact integro-differential evolution equation obeyed by the function J(p^,/i) is 



din n \_ jj,'^ J JO — 

(40) 

We encounter again the cusp anomalous dimension, and in addition a new function 7'', which 
has been calculated at two-loop order in ^31j. 

An important object in the derivation of the solution to this equation is the associated jet 
function j(L, /i), where L = ln(Q^/yU^). This function has originally been defined in terms of 
an integral over the jet function followed by a certain replacement rule 1221 • More elegantly, 
the associated jet function can be obtained from J by the Laplace transformation 



The inverse transformation is 

' c—ioo 



1 pc-\~ioo T 
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where the contour must be chosen to stay to the right of all discontinuities (i.e., c > 0). Using 
the evolution equation (^01) for the jet function, we find that the associated jet function obeys 
the RG equation 



(iln/i ' 



(in— ,/x 



2r 



cusp 



(43) 



which is local in and analogous to the evolution equation (^^1) for the hard function. The 
solution to this equation reads 



j(ln— = exp -4S'(/ii,/i) + 2ay(/ii,/i) 



2ar(Mi./i) ^2 



(44) 



where a^j is defined in analogy with (jg8|) . Given this solution one can readily derive the 
solution to the complicated evolution equation ()40|) for the original jet function by using the 
inverse transformation (|42p . The result is 



J(p^,/i) = exp — 4S'(;Uj,/i) + 2a^j{ni, n) 



-IeV 



m Jo 



dp'- 



p 



,/2^1-^7 ' 



(45) 



where 7] = 2ar(yUi,/i). This solution is valid as long as 1] > 0, which implies that fi < fii. 
Equation (jl^j) is analogous to the solution for the evolution equation of the 5-meson shape 
function found in |211E21 using a technique developed in pS] . 

Using the connection between J and j implied by Laplace transformation, it is possible to 
derive an even more elegant expression for the jet function J(p^,/i), which does not involve 
an integral and which is valid for both /i > /Xj and fj, < fii. The result relates J to the 
associated jet function j evaluated at the scale yUj, where it can be computed using fixed-order 
perturbation theory. We obtain 



J(p^/i) = exp -4:S{fii,fi) + 2a^j{fii,fi) j{d^,fii 



1 E_ 



2\V 



-lEV 



(46) 



where 5^ denotes a derivative with respect to the quantity t]. The star distribution is defined 



as 



dp' 



P^ lyU^ 



2 \ ^' 



f{p' 



+ 



/(O) (Q 



2\ ^ 



p. 



(47) 



where /(p^) is a smooth test function. The subtraction term involving /(O) is required only if 
7] < 0. For small rj, the above definition implies the expansion 



1 l_ 

p2 I ^2 



2 \ ^' 



6{p' 



+ 



+ r] 



1 f 



(4J 



The singularity for ^ is removed by the factor l/r(?7) in (j46j) . In the form given above, the 
expression for J{jP',n) holds as long as 77 > —1, which is sufficient for all practical purposes. 
For even smaller values of 77 it would be necessary to perform further subtractions in ()47|) by 
using the double-star distributions introduced in 
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3.3 Matching conditions and anomalous dimensions 

To evaluate the resummed hard and jet functions at a common factorization scale /i requires 
perturbative expressions for the matching functions Cv{Q'^,fJ'h) and j{L,fii). We extract the 
hard coefficient at a scale /i/t ~ Q in the first matching step, and the associated jet function 
at a scale /Xj ~ QV^ ~ ^ in the second. In this way, the matching functions are free of large 
logarithms and can be reliably computed in fixed-order perturbation theory. We also need 
perturbative expressions for the anomalous dimensions Fcusp, 7^, and 7"^. 

The hard matching coefficient Cy((5^,/i) is extracted in the first matching step, when the 
vector current in full QCD is matched onto an effective current built out of operators in SCET. 
To obtain an expression for the Wilson coefficient one must compute, at a given order in as, 
perturbative expressions for the photon vertex function in the two theories. The calculation 
is simplified greatly by performing these calculations on-shell, in which case all loop graphs 
in the effective theory are scaleless and hence vanish. The bare on-shell vertex function in 
QCD (called the on-shell quark form factor) has been studied extensively in the literature. 
The form factor is infrared divergent and can be regularized using dimensional regularization. 
The bare form factor at two-loop order was calculated long ago jlH |12l HSl HI] , and recently 
the infrared divergent contributions have even been computed at three- loop order jlHj. When 
the (vanishing) SCET graphs are subtracted from the QCD result, the infrared poles in 1/e 
get transformed into ultraviolet poles. To obtain the matching coefficient we introduce a 
renormalization factor Zy, which absorbs these poles. We then compute 

CviQ', fi) = lim Zy\e, Q\ fi) F^^U^, Q^) , (49) 

where on the right-hand side we must also eliminate the bare coupling constant in favor of the 
renormalized coupling as{fi)- At two-loop order, we find (with L = \n{Q^/ fx"^) and = as(Ai)) 

CviQ^f^) = 1 + ^ + 3^ - 8 + y) + Cf (^)' [CfHp + CaHa + TpUfRj] , (50) 

where 

^. = =r-3^^+ I ^-^1 1^^+1-^-^ + 24(3)^ + ^ + 1^-^-30(3, 




/ 45 
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'2545 
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Ha = —L'+\ -— + — +—- + —-- 26C3 U 



_ 51157 3377r2 Utt^ 313 
~ ~648 108~ ^ ^5~ ^ ~9~ 

4 . 38 . / 418 47r2\ 4085 237r2 4^ 

Hf = — L^H L^+\ LH \ h-Cs- (51) 

^ 9 9 V 27 9 / 162 27 9 ^ ^ ' 



This result agrees with the corresponding expression given in [21]. The anomalous dimension 
of the vector current in SCET is obtained from the coefficient Zy"* of the 1/e pole term via 
the relation ^ 

rcusp(«.) In ^ + i^'ias) = 2as ^ 4'^ iQ\ f^) ■ (52) 
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Using the results of [IHI the anomalous dimension can be extracted at three-loop order. We 
reproduce the well-known expression for the three-loop cusp anomalous dimension Fcusp |46j . 
For the quantity 7^, we obtain 



7^' (a,) = i - (4.68 - 0.95n/) — - (23.43 - 4.05n/ + 0.029n1) — + . . . . (53) 

The exact analytic expressions for the expansion coefficients are given in the Appendix. 

The two-loop expression for the jet function has recently been obtained in ^1] starting 
from expression (j2(Jj) . by which the jet function is expressed in terms of a two-point vacuum 
correlator in full QCD. Using those results, the two-loop matching condition for the associated 
jet function is found to be^ 



CpOis f^r2 or , ^ 271^ \ , ^ /O^N 



j{L, /x) = 1 + ^ i^2L'' -3L + 7-^j+CF i^^f [CfJf + CaJa + TpUfJf] , (54) 
where 

^ 9 V 18 3 J V 54 9 ^7 

53 1 29 1 557r2 377r^ 

+ I8C3 , 

648 36 180 ^ ' 

j,.?i3_55i^+fl5i_i£!^i_125I + l^. (55) 

^ 9 9 V 27 9 ; 162 9 ^ ^ 

The anomalous dimension kernel entering (pn|) has been calculated at two- loop order [3T]. In 
Section 13.51 below, we will show that the difference (7'' — 7^) multiplies the S{1 — z) term 
in the Altarelli-Parisi splitting function Pq«_g(z), which has recently been calculated at three- 
loop order jlB] ■ When combined with this relation serves as a cross-check of the two-loop 
result obtained from the direct calculation in [31], and further it can be employed to extract 
the three-loop coefficient of the jet-function anomalous dimension. We find 
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7^(a,) = -— - (0.364- 0.556n/) j - (3.18 - 1.33/2/ + O.OllnJ) (^—j +.... (56) 
The exact analytic expressions for the expansion coefficients are given in the Appendix. 



3.4 Resummation of large logarithms 

We are now ready to write down a resummed expression for the structure function F2'^{x^ Q^), 
valid to all orders in perturbation theory and at leading power in {1 — x) and K^q-^I^x- The 

^Our function j[L, jj) should coincide with the object M. n,dis derived in j5J after the substitution L — L. 
Comparing the two expressions, we disagree with the signs of the two- loop 0{L) terms with color structures 
CfCa and CpTpUf, and with the two-loop constant terms with color structures Cp and CfCa- 
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result is 



2 \ -2ar(Mh,iU/) 



X exp 



AS{nh,l^f) - 4:S{iJ,i,iif) - 2a^v{nh,l^f) + 2ay(/Xi,/i/) 



p-lEV 



r{r)) 



^(e,/^/), (57) 



where rj — 2ar(//i, as above. To leading power, we could approximate {^/x — 1) — > 
— x), but we prefer to keep the full x dependence in our numerical studies below. The 
"factorization scale" fif = is, by definition, the scale at which the parton distribution 
function is renormalized. 

The Sudakov exponent can be simplified using the general relations 



S(ijLi, 112) + 'S'(//2, A^s) = 'S'(//i, /ig) + In — ar(//2, A^s) ■ 

A'2 



Introducing the short-hand notations 

we find after a straightforward calculation 



Qj^tj) QjryJ QjryV ^ 



(58) 



(59) 



FT{x,Q')^Y.^l\Cy{Q^^xu)\' 



-2aT{lJ.h,IJ-i) 



exp [4S'(/Xfe, /Xj) - 2a^v{i^h, l^i) 



X exp 



2a^,A(/Xj,/i/)J j(^ln — + 

A^i 



1— ?7 



(60) 



The remaining integral can be performed by noting that, on very general grounds, the 
behavior of the parton distribution function near the endpoint can be parameterized as 



(61) 



where 6(/x/) > 0. We will see in the next subsection that this functional form is preserved 
under RG evolution. Defining a K factor by the ratio 

Fr(a;,g2) 



K{x,Q\i,f)^ 



J2gelx(f)f{x,iif) ' 



(62) 



we now obtain our final expression 

K{x,Q^^if)^\Cv{Q',^iH)\'(^^' 



-2ar {^ih,^^i) 



exp 



AS{nh,lJii) - 2a^v{iJih,iJii) 



(63) 



X exp 



j( In — 

V X 



X 



g2 1 - a; ^ _ X e 



r(i + %/)) 



r(i + %^) + 7y) ' 
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Table 1: Different approximation schemes for the evaluation of the resummed 
factorization formula ()63|) 



RG-impr. PT 


Log. Approx. 


Accuracy ~ a^L^ 


r 

cusp 


7\ Y 


Cv,J 




LL 


n + l<k<2n (a^^) 


1-loop 


tree-level 


tree-level 


LO 


NIL 


n<k<2n (a°) 


2-loop 


1-loop 


tree-level 


NLO 


NNLL 


n — 1 < k <2n (a^) 


3-loop 


2-loop 


1-loop 


NNLO 


NNNLL 


n-2<k<2n {a^) 


4-loop 


3-loop 


2-loop 



where as before r] = 2ar(/ii, /i/). In this expression, the dependence on the two physical scales 
and M"^ = (neglecting the nucleon mass) is completely explicit. Our exact result is 

independent of the scales fit and /Xj, at which QCD is matched onto the intermediate and final 
effective theories, SCET(/ic, c, sc) and SCET(c, sc), respectively. In practice, a residual scale 
dependence remains once we truncate the perturbative expansions of the various objects in the 
factorization formula. The final answer simplifies further if we choose the "natural" values of 
the two matching scales, = Q and /Xj = Mx- However, we prefer to vary the matching scales 
over some reasonable range and take the variation of the results as an estimate of higher-order 
perturbative effects. Note that by definition the K factor does depend on the choice of the 
factorization scale fif. It is conventional in the literature on DIS to identify fif with the hard 
scale Q. However, from the point of view of an effective field theory it would be more natural 
to choose fif below the intermediate matching scale fii ~ QV^ "~ A typical choice would be 
/X/ of order a few GeV, independent of the dynamical variables x and Q. 

In (|Um) and (jU^ we have accomplished the complete resummation of threshold logarithms 
for F2 directly in momentum space. That the final answer is a convolution (rather than a 
product) of a hard-scattering kernel with the parton distribution function is reflected in the 
non-trivial dependence on the hadronic parameter 6(/i/) describing the large-^ behavior of 
0g'^(^,/i/). Our factorized expression for the DIS structure function is very similar to that 
for the B Xg'j decay rate derived in [23]. Although the hard functions and soft matrix 
elements differ, the jet function is the same in both cases. An important advantage of our 
momentum-space approach is that in the limit where the two matching scales are set equal 
to the factorization scale, /i/i = yUj = /i/, the resummed results (jUUI) and (jUHjl automatically 
reduce to the corresponding expressions valid in fixed-order perturbation theory (expanded 
about X = 1). Consequently, it is straightforward to match our resummed expressions onto 
fixed-order calculations valid outside the threshold region. 

The right-hand sides of (j6(Jj) and (j63|) can be evaluated at any desired order in resummed 
perturbation theory. Table Q shows what is required to obtain different levels of accuracy in 
the perturbative evaluation of the result. In this work we adopt the counting scheme of RG- 
improved perturbation theory, where at leading-order (LO) one includes all terms of order 1, 
at next-to-leading order (NLO) one includes all terms of order a^, and at next-to-next-to- 
leading order (NNLO) one includes all terms of order a^. We count the large logarithms 
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Figure 5: Comparison between fixed-order (dashed) and resummed results (solid) for the K 
factor. The green curves are the LO result, red NLO, black NNLO. For the resummed result, 
we set iih = Q, fJ^i = Mx, fij = Q, and b{fif) = 4. The fixed-order result is obtained by setting 
all scales equal to /i/. 

L e {In fih/f^i, ^^fJ'i/f^f, ln(l — x)} like l/ag- In the literature on DIS, the LO approximation 
is also referred to as the next-to-leading logarithmic (NLL) approximation, the NLO result 
is referred to as the next-to-next-to-leading logarithmic (NNLL) approximation, etc. The 
leading logarithmic (LL) approximation is listed only for completeness, as it neglects terms 
that are parametrically much larger than 1. 

In Figure we compare the fixed-order calculation of the K factor with the resummed 
result for Q = 5 GeV and Q = 30 GeV. For the resummed result we use the default choice of 
scales fih = Q, f^i = Mx = Q\/^^ and take the asymptotic form of the parton distribution 
(16 ip with b{fif) = 4 in both cases. Following common practice we choose fif = Q for the 
factorization scale. In this case the quantity rj < 0, and because of the factor (i^)*? in 
(I63|) the resummed results diverge as x approaches 1. The figure illustrates that higher-order 
corrections become important as a; — *• 1, and that fixed-order perturbation theory is no longer 
adequate in this limit. The magnitude of the K factor can be reduced by adopting a lower 
choice for the factorization scale, which is more in line with the philosophy of an effective 
field-theory approach. For example, we may consider taking fif ^ Mx{x = 0.9) ~ <5/3, 
corresponding to a typical hadronic invariant mass in the endpoint region. The corresponding 
results are shown in Figure El We observe that with such a choice of the factorization scale 
the K factor takes more moderate values, and also that the effects of resummation are less 
significant. 

In Figure [71 we show the scale dependence of the result obtained by varying the hard and 
intermediate scales by a factor of 2 about their default values. The figure shows a dramatic 
reduction in scale uncertainty when going from LO to NNLO. It also suggests that varying 
the two matching scales individually by a factor of 2 may overestimate the perturbative un- 
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Figure 6: Same as Figure El but with a lower choice of the factorization scale. Specifically, 
we take /i/ = 1.5 GeV for Q = 5 GeV (left), and /i/ = 10 GeV for Q = 30 GeV (right). 




Figure 7: Scale variation of the K factor at Q = 30 GeV. The light-gray band is obtained 
by varying Mx/2 < fii < 2Mx, while the dark-gray band arises from varying the hard scale 
Q/2< i2h< 2Q. We set /i/ = 30 GeV and 6(/i/) = 4. 
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Figure 8: Same as FigurelHl but now including the suppression factor (1 — xY with = 4 
from the parton distribution function. 



certainty, because the higher-order results lie near the center of the large band obtained by 
varying the renormalization scales in the low-order ones. A variation of the scales by a factor 
of a/2 better represents the uncertainty in the present case. Furthermore, it seems reasonable 
to perform the scale variations in such a way that the hierarchy of scales /i^ > fii is preserved. 

We stress that the resummation of large logarithms accomplished in fl63p is under pertur- 
bative control as long as (1 — x) ^ AQ^jy/Q"^, since only then is the intermediate matching 
scale jji ~ gVl — X ^ Aqcd a short-distance scale. Physically, this condition is equivalent to 
saying that the final-state jet can be treated in an inclusive way using a partonic language. 
Numerically, we can assume that perturbation theory at the jet scale breaks down in the region 
where Mx < 1 GeV. We illustrate this boundary in our x-space results with the gray band 
in Figure El For Q = 5 GeV, our approach is valid as long as x < 0.96. For Q = 30 GeV, it 
extends all the way to x < 0.999, so that in this case the band is not visible on the scale of the 
plot. While our theoretical description breaks down very close to the endpoint, we note that 
weighted integrals of the jet function over an interval x G [xq, 1] can be calculated starting 
from (jUUI) as long as xq is in the short-distance domain. 

Experimentally, structure functions at large x are very difficult to measure, because of the 
rapid decrease of the parton distribution function as x — > 1. This is illustrated in Figure |Hl 
where we include the suppression factor from the parton distribution. Because of this strong 
suppression, there are no measurements of the non-singlet structure function for x > 0.9 
|471 1^ H^. The experiments that probed the highest x- values were fixed-target experiments 
in the 1970s and 1980s at SLAG [50] and the BGDMS experiment at GERN [H]- Newer 
experiments only cover x < 0.65. As a consequence, the threshold resummation in DIS is 
currently of limited phenomenological importance. However, the resummation is relevant for 
W or Higgs production at hadron colliders, which can also be analyzed with the methods 
developed here. 
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3.5 Parton evolution near the endpoint 



The easiest way to derive the evolution equation for the parton distribution function in the 
limit a; —> 1 is to use that the factorized expression (j^Uj) for the structure function F2 must be 
independent of the arbitrary renormalization scale /i, and to combine this information with 
the known scale dependences of the hard and jet functions, given in (jHBjl and (jlUI) . This yields 

d /•! ^^''i^' a) 

= 27na.)0-(e,/i) + 2r,.,p(a.) 

= y^'f ntrH^)0r(|/^)' (64) 

where 

p(endpt)(^) = ?^f^ + 27^(a.) 5(1 - z) (65) 

is the z ^ 1 limit of the Altarelli-Parisi splitting function Pq^g{z), which is known from 
direct calculation at three-loop order The asymptotic form of the splitting function 

near the endpoint given above holds to all orders in perturbation theory, up to corrections of 
order {1 — z). Recall that the anomalous dimension 7*^ was defined as the difference of the 
anomalous dimensions 7"' and 7^ of the jet function and SCET current, see (jSHI)- Relation 
(jU^ thus provides a check of our two-loop results for these anomalous dimensions, and it 
furthermore allows us to deduce the value of the three-loop coefficient 72 given in relation 
(jnZ|) of the Appendix. 

The exact solution to the evolution equation ()64|) can be found in analogy with ()45|) . It 
reads ^ 

/^/) = [2 V (/^/' /^o)] di' _ , (66) 

where this time a = 2ar(yUj, /xq)) and fiQ denotes the scale at which the initial condition for 
(f)"^ is given. For the hadronic parameters A/" and b governing the asymptotic behavior of the 
parton distribution function in this relation implies 



Kf^f) = KfJ'o) + 2ar(/i/, /io) , 

Afifif) = Ar(/io) exp [2a,,(/i/,/io)J ^,,H,f) ^ ^^^^^^^ ■ (67) 

These evolution equations ensure that the fif dependence on the two sides of relation (j63|l is 
indeed the same. The first result is particularly simple and interesting. Since ar(/i/,yUo) > 
for fif > yUo; it follows that the coefficient b{n) increases with fi, a fact incompatible with the 
naive counting rule result 6 = 3 In other words, such a counting rule could possibly 

hold only at a specific renormalization point. 



4 Connection with the standard approach 

The conventional approach to threshold resummation in DIS proceeds via moment space and 
inverse Mellin transformations mi2]- The purpose of this section is twofold; first, to show that 
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our momentum-space resummation is formally equivalent to the conventional resummation 
order by order in perturbation theory, and second, to point out and quantify the theoretical 
and numerical differences that appear in applications to physical quantities such as the DIS 
structure function. 

We first recall some details of the conventional approach to threshold resummation in 
moment space. By taking moments of the structure function ^2"*^, convolution integrals such 
as (jllj) or can be brought into product form. The traditional way of writing the result is 



7-1 ns 



dxx 



N-l T^ns 



(6^ 



where the moments of C{Q'^, z, fi) and (f)q^{C,, /i) are defined in analogy to those of F2^{x, Q^). 
For large N, the function Cat is then split into two pieces according to 



Cjv(g',/i/) = (7o(Q',/U/) exp [Gjv(Q',/i/)l +^(]^ 



(69) 



where Qq contains all A^-independent contributions, while the function Gn contains logarithms 
of the form In'^ N. The limit a; — > 1 in momentum space corresponds to the limit — »• cxd in 
moment space, so this formula achieves the exponentiation of large threshold logarithms. The 
resummation exponent Gat is written as 



GNiQ ,IJ'f)= / dz 



1 



Ag{a,{k)) + B,(^as{QVT 



(70) 



where the functions Ag and Bg are universal radiation factors determined by matching with 
results from fixed-order perturbation theory. 

We shall now derive an equation relating the objects go, Ag, and Bg in (jU^ and (fTUj) to the 
matching coefficients and anomalous dimensions defined in effective field theory. We begin by 
transforming the factorization formula into Mellin space, obtaining the product form 



F2:NiQ') = |Cv(Q^/i/)|v^(g^/i/) ^ e: 



?Civ+i(/^/) 



(71) 



which is valid up to corrections in The Mellin-transformed jet function is defined as 



7V-1 



J^(g^^) = dp'h-^] Jip^f,). 



It was shown in that for large the jet-function moments Jn are given by 



(72) 



(73) 



and hence obey the same evolution equation pHjl as the associated jet function. Using this 
connection along with the results derived in Section IHl the resummed coefficient function Cn 
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in (jHS|) can be written as 



exp 



4S{iJ,h,Hi) - 2a^v{fih, fJ-i 



X exp 



~ 1 
2a^^(/ii,/i/) - 21nA^ar(/ii,/i/)] j(ln-^,/ii) + • (74) 



We now adopt the "natural" scale choices fJh = Q and /Xj = Q/vN, which are implicit 
in most treatments of threshold resummation in the literature. This allows us to compare 
with the standard expression (jEHI), but as we will discuss at the end of this section, this scale 
choice becomes problematic when the expressions for the moments are transformed back to 
a;-space. Next, we express the RG functions S'(/ii,/i2) and a„(/ii,/i2) defined in in terms 
of integrals over the appropriate anomalous dimensions. After a straightforward calculation, 
this leads to 



^o''''^(Q',/i/) = \Cv{Q\Q)\'' J(0,Q) exp 



-\nN 



1 in / /, 7/ /, d\nj(0, k) 

In r,,,p(a.(/.)) - r'Mk)) - -j^^ 



1 



— ^cnsp{0is{k)) 



(75) 



where we have defined the split between the two terms such that the expression for Gjsi 
obtained in the large- limit vanishes for iV — > 1. 

Our next task is to bring the exponent G^r from the standard result ()7Up into a form 
resembling the SCET result fl75|) . Since the running coupling as{k) depends on its argument 
logarithmically, a helpful identity is (for integer A; > 0) 



1 ^JV-l _ 1 



dz ■ 



1 - z 



In^il - z) 



^ 4+i(lni)+of^ 
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where 



In{x)=d: [e^("+^^)r(l + e) 



(76) 



(77) 



are n-th order polynomials defined in 
large N 



With the help of these relations we find that for 
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Ag{asik)) - Bgia^ik)) 
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InN 
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where 



AG(/z) = E [Af-^'Mi^)) + 5f («.(/i)) 

k=i y'^ + 



(7J 



(79) 
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and denote the n-th derivatives of Aq and Bq with respect to In jj?. The perturbative 

expansion of AG starts at order a^. Contrary to the SCET expression in (f7H|). the result (f7S|) 
does not vanish for iV — 1. The overall normalization of Gi\f{Q'^, ^f) is a matter of convention, 
since it can be absorbed into go{Q'^,fif)- Taking the difference in normalization into account, 
the two definitions underlying (fTSj) and (fTHj) are connected by 



^o(Q^/^/) 



^SCET 



(Q^fij) exp [-AG(Q)] 



^0) 



At the expense of a proliferation of 7^; terms in the perturbative expressions, one can equally 
well normalize = 0. This normalization condition is adopted, e.g., in j5. 

Equation (f7S|) is consistent with (f73|) if we identify Ag{as) = rcusp(«s) with the cusp 
anomalous dimension, and furthermore require that 



Bq^asifx)) + 



din 

This formula can be rearranged to read 



d\iaj{0,fi) 
din fi'^ 



r(v) - = 



cusp 



^1) 



>.) , (82) 



where 



d 



and the differential operators are defined by their Taylor expansions, 
perturbation theory we obtain^ 



(83) 

Evaluating (|82|l in 
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2C3ro7n + — 



(84) 



where the one- and two- loop matching coefficients B'q^ and Bq'^ have been calculated in [3T]. 
Computing the first three Bg^n coefficients using the three-loop result for the anomalous di- 
mension 7"^ given in the Appendix, we find agreement with the expressions derived in 
Ref. identified the function Bg with the jet-function anomalous dimension 7"^, which is 
incorrect already at two- loop order. 

Equations (jH^ and (j5^ provide the desired relations between the function Bg of the 
standard approach and the field-theoretical objects defined in the effective theory. Obviously, 
the connection between the various objects is highly non-trivial. This explains, perhaps, why it 
has proven difficult in the past to translate between the standard formalism and the approach 

•^A relation similar to H82(l has been derived in 21 ; however, there is a typo in the last equation in (75) of 
that paper, which is the analog of our relation between B^ ^ and 7^. 
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based on SCET. The deeper reason is that in the conventional approach the RG evolution 
equations of SCET are replaced by a different set of partial differential equations P3I21, whose 
solution is equivalent to our solution but not structurally identical to it. In particular, there 
are theoretical and numerical differences between the moment-space and momentum-space 
resummation procedures. Some of these are explicit in the particular form of the resummation 
exponent ()70|) obtained in moment space, and some become apparent only when performing 
the inverse Mellin transform. We conclude this section by examining these differences in more 
detail. 

A troublesome feature of the conventional moment-space approach is that the integrals 
over the coupling constant in the resummation exponent run over the region where ^^(/i) is 
evaluated at very small values of /i. To leading order, the coupling behaves as 

An 

and becomes infinite at the scale n = Aqcd- When the integration variable z in (f70|l ap- 
proaches 1, the resummation exponent becomes sensitive to this Landau-pole singularity in 
the running coupling. As a result the integral is ambiguous, since one can arbitrarily choose 
a prescription for dealing with the pole. We can estimate the magnitude of the ambiguity by 
taking the difference of the z-integral evaluated above or below the Landau pole in ()70|) . The 
result is 

AG„ = -|^(r„ + -,„^)(iV-l)^ + o(^!^), (86) 

which is of the form of a power correction of order Aqqj-,/M|^. Note that this ambiguity 
never appears in the momentum-space formulation, and should therefore be interpreted as 
an artifact of resummation in moment space. As stressed earlier, the Landau-pole ambiguity 
does not imply an infrared renormalon-pole ambiguity of the same strength Aqqj-,/M|^. To 
show that Gat is indeed affected by a corresponding renormalon pole, one needs to evaluate 
the exponent in the large-/3o limit, a fixed-order truncation of this quantity is not sufficient 
[7j. On general grounds, one expects anomalous dimensions to be free of infrared renormalons, 
so that the renormalon poles enter only through the associated jet function j{0,fi) in ()82j) . 
In the effective theory, renormalons affect only the matching coefficients, Cy and j, and will 
always be commensurate with power-suppressed operators. RG evolution, on the other hand, 
is driven by anomalous dimensions which are free of renormalons. 

Further differences become apparent when studying the inverse Mellin transformation 
needed to obtain the physical momentum-space results from the moment-space expressions. 
While our result (|63|) obtained directly in x-space is completely analytical, the inverse Mellin 
transform can only be performed numerically, by evaluating the integral 

iTTl Jc-ioo 

The Mellin inversion is actually ambiguous, because the expression for Cn has a Landau pole 
for large N. We deal with this pole by adopting the so-called minimal prescription [5], which 
amounts to excluding the Landau pole from the integration contour by choosing the constant 
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Figure 9: Comparison between Mellin-inverted moment space results (dashed) and results 
obtained in directly in x-space (solid). The green curves are the LO result, red NLO. The 
black lines are NNLO results and are visually indistinguishable from the NLO curves for 
Q = 30 GeV. We set /x/t = /i/ = Q, and &(/i/) = 4. For the intermediate scale, we choose 
/ij = Mx in momentum space and fii = Q/ y/N in moment space. 



c smaller than the value of N at which the pole occurs. Even with this prescription, the 
numerical integral is not well behaved in the limit x —>■ 1, since the damping of the integrand 
becomes weaker and weaker as x approaches the endpoint. In Figure IHl we compare the results 
for the x-space structure function obtained through numerical Mellin inversion with those 
obtained directly in momentum space ()63p. One source of numerical differences arises because 
the relation (fTSj) is only approximate,^ so that the solution to the RG equation for JAr(Q^,/i) 
receives corrections which are suppressed as while our momentum-space solution (j46|) 

is exact. Another is that the default choice of the intermediate scale fii is different in the 
two approaches. The numerical differences are noticeable for smaller values of Q, but become 
negligible at Q = 30 GeV. 

In the effective-theory result for the moments, the Landau pole in the inverse Mellin trans- 
formation can be avoided by performing the inversion to x-space with the appropriate scale 
choice for momentum space, ~ QV^ — x, instead of /ij = Q/Vn. The freedom to choose 
the scales as appropriate for the quantity under consideration is an important advantage of 
our approach. The Landau-pole ambiguity in the Mellin inversion is not the only problem that 
arises from the fact that the scales cannot be varied in the standard resummation formalism. 
An additional difficulty was pointed out in [5 . To illustrate it, let us consider the structure 
function at the leading logarithmic level, even though this is not a consistent approximation 
in RG-improved perturbation theory. Our result then reduces to 

K{x,Q'^,fif) = exp [4:S{fih,fJ'i) + 2ar(Aii,/U/) ln(l - x)] , (88) 
^The exact form of the RG equation obeyed by the jet-function moments can be found in pi] . 
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where we have approximated ^— ^ !^ 1 — x. From (|63|) . we see that we have to choose fit ^ Q 
to make the double logarithms in the perturbative expansion of the hard matching coefficient 
Cy(Q^,/i) small. Similarly, to avoid the appearance of large logarithms in the associated jet 



function j(\n—f,fi), the choice fii ~ Mx is mandatory. Let us now look at the structure 



function integrated over the endpoint region 

m^,Q')= C dyFr{y,Q'). (89) 

Jl-x 

In this case, the appropriate choice of the intermediate scale for integral J-'2^{x,Q'^) is fii ~ 
Q\/l — X, as can be checked by explicitly performing the integral over ()63|) . If one instead 
chooses the scale /ij to avoid logarithms on the level of the integrand, then the integral 



becomes singular. To see the problem, we set fif = fih = Q, /^i ~ QV^ ~ V ^^id, for illustration 
purposes, approximate the Sudakov factor by expanding it to leading order around fixed 
coupling as(Q), as was done in [5]. The integral (|H^ then becomes 



2 2 

-aln2^ + 2aln^ \nil-y) 



^-(x, Q') = f dyY^e^y 0^^(i/, Q) exp 

Jl-X g 

= C rf2/Ee'2/C(2/'/^/)exp[aln'(l-y)l, (90) 

Jl-X g 

with a = Fq ° g^^^ . Because the exponential factor grows faster than any power as ?/ — > 1, this 
integral diverges. Its expansion in a is an asymptotic series with factorially growing terms. 
As was shown in ^ the ambiguity associated with the non-integrable singularity for ?/ ^ 1 is 
of order 

^_l/4a _ /Aqcd\4§^ _ /AqCD\1-4 ^g^^ 



Q ^ ^ Q ^ 

for = 5. In [S^ it was shown that the above divergence does not occur if the Sudakov resum- 
mation is performed in moment space and the inverse transformation is performed exactly, 
without dropping subleading logarithms ln(l — x). From this, the authors concluded that 
the appropriate place to perform resummations is moment space and that leading logarithmic 
resummations in x-space are problematic. Our analysis shows that it is simply a bad choice 
of scale that produces the problem of the spurious power correction: the usual moment-space 
formalism produces logarithms In^ in the Sudakov exponent, which translates into ln^(l — x) 
at leading logarithmic accuracy, which in turn causes the problem in ()90|) . However, the proper 
way to perform the calculation is to keep the matching scales arbitrary and choose them such 
that the final result of a given calculation does not contain large logarithms. This avoids 
the above problem as well as the occurrence of Landau-pole ambiguities in inverse Mellin 
transforms. 

We hope that the above discussion helps to overcome the misconception that Mellin mo- 
ment space is the "correct place" to perform the threshold resummation, and that resummation 
in x-space leads to inconsistent results. Quite to the contrary, the final analytical formulae we 
obtain in momentum space are simpler than those derived in moment space, they are free of 
spurious, unphysical power ambiguities and, as Figure El shows, the perturbative expansion in 
x-space exhibits a better apparent convergence. 
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5 Summary and conclusions 



We analyzed DIS in the threshold region x — > 1 using SCET. With a detailed analysis in the 
effective theory, we rederived the standard QCD factorization theorem ^1 121 El for the non- 
singlet structure function F2^{x, Q^). While this process had been investigated in the effective 
theory before [HI UHl 1111 12II1 1211 122|, we argued that previous studies were incomplete. Our 
analysis resolves the issues left open in these papers. We agree with that in a diagrammatic 
analysis momentum modes with low virtuality appear. Their presence is a consequence of the 
fact that in the limit x — > 1 one parton carries nearly all the momentum of the nucleon 
and the characteristic scale associated with the target remnants is m^(l — x), at least in a 
perturbative analysis. However, here we have shown that these modes do not translate into 
non-perturbative Q dependence in the parton distribution function. In ^7] it was argued that 
such momentum regions do not contribute to the effective-theory calculation, and in [20] that 
they would be screened away non-perturbatively. Since these target remnants are part of the 
endpoint parton distribution function, it is incorrect to exclude them from the effective-theory 
factorization analysis, and there is no need to invoke a mechanism which eliminates them. 
For the same reason, there are no extra "soft" contributions outside the parton distribution 
function, contrary to what was postulated in [TH] . 

With the factorization theorem at hand, we then performed the resummation of large log- 
arithms by solving the RG equations of the effective theory. Our result involves three scales: a 
hard matching scale fih, a jet scale /ij, and the factorization scale /i/. By choosing the pertur- 
bative scales to satisfy /ih ~ Q and /Xj ~ Mx, we avoid the presence of large logarithms. This 
approach has several advantages compared to the standard resummation technique for DIS. It 
enables us to derive a simple analytic expression for the resummed structure function directly 
in x-space, thus circumventing the problems associated with moment-space resummation. We 
can also estimate the higher-order perturbative uncertainties by varying the matching scales. 
It is trivial to recover the fixed-order result by setting all scales equal, fit = fJ^i = fJ^f. This 
makes it straightforward to combine our results with fixed-order calculations valid away from 
the threshold region. 

An advantage of the effective-theory approach is that the resummed results are free of 
Landau-pole ambiguities. In the standard approach, these appear twice: in the resummed ex- 
ponent in moment space and also in the Mellin inversion back to momentum space. Since we 
perform the resummation in momentum space by integrating out the higher scales and using 
RG evolution to go to lower scales, our expressions do not involve the strong coupling constant 
evaluated at scales below the minimum of the intermediate matching scale and the factoriza- 
tion scale at which the parton distribution function is renormalized. Therefore, Landau-pole 
ambiguities do not arise at any finite order in perturbation theory. We showed that our re- 
sults are formally equivalent to the standard ones order by order in perturbation theory. This 
allowed us to relate the radiation function Bg to a combination of the anomalous dimension 
7^^ of the jet function and effective-theory matching coefficients. The two objects Bq and j"^ 
are identical at leading order, but beyond this the relation is highly non-trivial. 

Since the parton distribution functions fall off very rapidly near x — > 1, it is experimentally 
challenging to measure structure functions at large x. For this reason the amount of avail- 
able experimental information near threshold is very limited. However, because of its relative 
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simplicity and since the perturbative quantities are known with high precision, the threshold 
rcsummation for DIS has provided us with an ideal setup to develop our formalism. In the fu- 
ture, we plan to use the same approach to perform resummations in other, phenomenologically 
more relevant situations. 
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Appendix 



The exact solutions to the RG equations (jHHj) can be evaluated by expanding the anoma- 
lous dimensions and the QCD /3-function as perturbative series in the strong couphng. We 
work consistently at NNLO in RG-improved perturbation theory, keeping terms through order 
dig in the final expressions for the Sudakov exponent S and the functions ar, a 
We define the expansion coefficients as 
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and similarly for the other anomalous dimensions. In terms of these quantities, the function 
ar is given by 
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The result for the Sudakov factor S is more complicated. We obtain 
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where r = Q;s(/i)/as(z/). Whereas the three-loop anomalous dimensions and /5-function are 
required in (jHSI), the expression for S also involves the four-loop coefficients T3 and P3. 
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We now list expressions for the anomalous dimensions and the QCD /3-function, quoting 
all results in the MS renormalization scheme. For the convenience of the reader, we also give 
numerical results for n/ = 5. The expansion of the cusp anomalous dimension Fcusp to two- 
loop order was obtained some time ago |10j, while recently the three-loop coefficient has been 
obtained in [16J. For the four-loop coefficient F3, we use the Fade approximants derived in 
jl]. The results are 
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The anomalous dimension 7^ can be determined up to three-loop order from the partial three- 
loop expression for the on-shell quark form factor in QCD, which has recently been obtained 
in HSl . We find 
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The results for the expansion coefficients of the jet-function anomalous dimension are 
7„" = -3Cf = -4, 

,i ^ CI (-1 . 2.^ - 24C3) + (-^ - + 4OC3) + C.T,n, (f + 
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Finally, the expansion coefficients for the QCD /3-function to four-loop order are 
34 20 

(3i = —Cl-— CaTfUj - ACpTpUf ^ 38.6667 , (98) 
o o 

ft . 5551 cl . (2c^ - ?f - If ci) c. + f c.) r|„5 

^ 180.907, 

149753 /1078361 6508 . \ /50065 6472 . \ ^ 1093 , 



729 ^ 



149753 ^ /1078361 6508 , \ /50065 6472 ^ \ . 

« 4826.16, (99) 
where the value of is taken from [31] and corresponds to A^c = 3 and Tp = |. 
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